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Abstract 

We expose a functional integration method for the averaging of continuous products 
Pt of N x N random matrices. As an application, we compute exactly the statistics of 
the Lyapunov spectrum of Pt- This problem is relevant to the study of the statistical 
properties of various disordered physical systems, and specifically to the computation 
of the multipoint correlators of a passive scalar advected by a random velocity field. 
Apart from these applications, our method provides a general setting for computing 
statistical properties of linear evolutionary systems subjected to a white noise force 
field. 
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1 Introduction 



In this work we give a detailed exposition of a functional integral method for the averaging 
of time-ordered exponentials of N x N random matrices which has found several applications 
in the study of the statistical properties of disordered systems. 

The method was introduced by one of the authors in the N = 2 case in order to 
compute the partition function of the Heisenberg ferromagnet, and was thereafter applied to 
the study of one- dimensional Anderson localization || and to some problems of mesoscopic 
physics Q| . Later, the same technique |3|, [7| allowed to obtain analytical results in the 
problem of a passive scalar advected by a 2-dimensional random velocity field. The approach 
of refs. [EL [7| was then extended to the more general iV-dimensional case in ref . M . 

In the present work we expose the method in its full generality and show that it allows to 
compute exactly the statistics of the whole Lyapunov spectrum of the matrix P t describing 
the time evolution of a linear system subjected to a white noise force field. Such a spectrum 
is relevant || in the computation of the multipoint correlators of the passive scalar (as well 
as in the computation of the correlators of passive vectors and tensors); in this case the 
matrix P t describes the time evolution of particles in a turbulent fluid, linearized around a 
given trajectory. 

However, our setting presents a high degree of generality and therefore a wider range 
of applications. In particular, our formalism can be naturally applied to the study of TV- 
level quantum-mechanical systems affected by a random noise, and it is complementary to 



the supersymmetric approach [II, TE] to the problem of N/2 channels localization in a 
disordered conductor. 

2 Averages of time-ordered exponentials 

Let us start with the problem of computing gaussian averages of the form 

(nPt}) = jjj Z>*exp (S[X}) F[P t ] (1) 
where X(s), for < s < T, is a traceless N x N hermitian matrix, 

V ^= U YldX^dX^WdX^s) (2) 

0<s<Ti<j i 

is the Feynman-Kac measure, M is chosen in such a way that (1) = 1, 

S[X] = ±- F Tr X\s) ds (3) 



AD Jo 

and Pt is the time-ordered exponential 



P, T ■( f*X(s)ds) Mi 
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such that r(t) = iVo is the general solution of the linear problem r = Xr, r(0) = r and 

AP- 1 = X (5) 

We shall now introduce a set of "collective" integration variables in order to re-express (H) 
in a more tractable form. At the same time, we shall chose the new variables in such a way 
that (|3p still be quadratic and that the Jacobian determinant of the functional transformation 
be particularly simple. 

As a first step, let us Gauss-decompose the matrix P t : 

A=(l + 0(t))-exp(f(t))-(l + ^)) (6) 

where 



(7) 



<f>ij(t) 


= 0, i<j 


e tj (t) 


= o, i>j 




= Ti(t) Sij 




JV-1 


T N (t) 









Moreover, in order to ensure the equality Pq = 1 we shall impose 

0(0) = 0, 0(0) = (9) 

We would like now to re-express the "local" degrees of freedom Xij (t) in terms of the global 
ones <fiij(t),Tij(t),6ij(t). This can be accomplished by making use of the basis of the 
matrix algebra, which is defined by (ey)jy = SikSji and satisfies the commutation rules 

[eij, hi] = 5jk&u - SuCkj (10) 

In particular, one has 

e f eij e" f = e Tl " Tj % (11) 
From these relations the desired expression for X^ readily follows: 

Xij = <fiij + E ^ik&kj 
k 

+fiSij + (pyfj + fi4>ij + E 4>ikTk4>kj (12) 

k 

+ ^2(<j>ikAkj + A ik (j) k j) + <j>ikAki(f>ij 
k k,i 

where Ay = e Ti ~ T i J2k 0%k{hj + hi), 

4>ij = -<fiij + E (pik^kj - E 4>ik4>ki4>ij H — ( 13 ) 
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and a similar definition holds for (for any fixed N flT3|) is a finite sum, since is a nilpotent 
matrix; the same is true for 9). 

Substituting fll2|) in (|) one obtains 

1 - 1 N ( ■ N ~ ■ \ ( N ■ - ■ \ 

- Tr X 2 = - ^ r| + £ Uy + £ <fe<^ e^"- £ fy0 H + (14) 

z z i=i y \ fc=i / \fc=i / 

The form of (0) suggests the introduction of the new variables 

0y = E ^ ( 5 m + 0«) + <Pii eT] ~ Tl ) eT °~ rj > * < i ( 15 ) 

k,l 

<pij = 0, i > j 

so that 

1*1* 

Relation (|TjJ) can be inverted giving 

% = E ^^~ Ti (4* + <Pkie Tl - Tk ) (8n + e l3 ) X u (17) 
k,i 

where 

:<j as) 

Through (|H]) one can re-express the 6^ as functions of the new variables </>y, Tj and 0^ in a 
recursive way, thanks to the "triangular" form of the equation. For instance, for N = 3 one 
gets 

t 



Mt) = / 023(3)^ T1(SK2T2{S) ^, 
JO 

12 (*) = [ t A(s)ds, where A = (0 12 + ^s^K 2 ^ 1 , (19) 

JO 

= f[A( s )e 23 ( s )+fe(5) e - 2 ^ s )-^ s )]^ 



As a matter of fact, for any fixed % the N — i functions % can be expressed through the 
N — i — 1 functions j and the remaining variables by means of a single quadrature. This 
is an important point, since for practical calculations Pt has to be re-expressed in terms of 
the new variables <f), f, 0. 

We must now substitute X(s) as an integration variable in the functional integral ([[]) with 
the new variables 0(s), f(s), 0(s). Again using the commutation rules fllQf) , and renaming 
fj = Pi for convenience, we finally get 

k 

Xu = pi + ^2(4>ik4>ki - 4>ik4>ki), i = l,...,N (20) 
k 
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Xij = o,ji>j + pifaj + a,, (21) 
+ ^{4>ikPk4>kj + 4>ik4>kj + 4>ik<Pkj - 4>ik4>kj) 

k 

- J2(xjk4>ik4>ki4>ij + xii4>ik<j)ki<t>ij) 
k,i 

- Xmk4>ik4>kl4>lm4>mj, « > j 



k,l,m 

In ref. the AT = 3 case was explicitly considered. We observe that the matrix elements of 
X(t) = PtPf 1 transform as 

X tj {t) -> e^Xi^t) (22) 
under the global gauge transformation 

-> e^fajit), 4 (t) -> o /; (/). -> Pi (t) (23) 

where cr^ satisfies o"ifc + tx^- = cr^-, cr^- = — a^. 

As the last step, we must compute the Jacobian determinant of the functional transfor- 
mation (|2l|) . Notice first that through the shift given by 

TV 

4>ij — > <t>ij - Y 4>ik4>kj (24) 

k=l 

N 

Pi -> Pi - Yi^i^ki - 4>ik4>ki) 

k=l 

(which has Jacobian J7 7 = 1) one can reduce to the computation of 

Det ( (25) 
\S(f>SpS(f>J 

where X_ is the strictly lower triangular part of X. The Jacobian (p5| ) can be computed by 
means of the standard regularization procedure (see ref. 0) 

4>n = <f>(tn), Pn = p(t n ), <fi n = 0(*n) 

t n = hn, n = 1, . . . , M, ft = T/M -»• 0, M -> +oo (26) 

<w — i — , <Kt) -> — - — 



giving 



/tv-i a 



p 3 (s)ds) (27) 



Applying now the variable transformation X — > (0, p, 0) one sees that the functional inte- 
gral (D) reduces to 

Ml) = ^> / VtVpV}ex V (-S%p,i]yi(l + ty(l + 6)} (28) 



where 6 = 9[4>,p, 0] is obtained by solving (|17|) , f = J T p(s) (is, the /?i are constrained by 
jV' is the normalization factor and 



S ' = 2DJo UT,pl + T,^-^D^{N-k)p k \ds (29) 



1 f T 1 1 N N ~ 1 

k=l i,j k=l 

In ([!]) the functional integration is constrained to the surface 

r = {x ij (s) = x; i ( s ), o<s<t} (30) 

In refs. |2|, [7] it was shown, using the Cauchy theorem, that whenever T is holomorphic 
in the matrix elements X^ one can modify the integration surface To to the homotopically 
equivalent 

T 1 = { ( p ij (s) = ^ l (s),lmp i (s) = 0, 0<s<T} (31) 

without affecting the value of the integral. This means that in (|28|) may be regarded as 
the hermitian conjugate of 0. 

We would like to remark that expressions similar to (|2T|) were obtained in the framework 
of conformal field theory |16| . However, the explicit form of the variables and of the Jaco- 



bian J ', which are essential for any physical application of our method, were not computed 
there. 

3 The Lyapunov Spectrum 

We shall now define the Lyapunov exponents Xj, j — 1, . . . , N, according to the relation |§ 

Ai + ■ • • + \ k = ^ log Vol(P r vi, . • • , /W) (32) 

where v 1; . . . , is an orthonormal set of vectors generating a unitary k- volume. For the 
sake of definiteness we shall choose Vj = e,-, where {ej)i = 5ij. One has 



Vol(P T e u ...,P T e k ) 



\E« (33) 

\ a=l 



n 



where A a (Mk), a = 1, . . . , 4, h = ) are ^ ne ^ x ^ minors of the n x k matrix = 
[P T ei, . . . , P T e k }. Let <pj,Pj, j = 1,...,N be the vectors defined by {4>j)i = + 0^, 
( Pi )i = (Prk- Then 

i,k 
N 

= J^^ihj + dkj)^ (34) 
fc=i 

= eJ0j + 2 kj e Th (j) k 
k<j 
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From (]34]) and the multilinearity of determinants it follows 



A a {M k ) = A a [ Pl ,...,pJ 

= A a [e n e 7 *^] (35) 



where Ai is the minor obtained from the first k rows of M k , are strictly positive integers, 
and elf are integer coefficients with c\' k = 0. One has then 



\ 1 + ... + \ k = I[ Tl (T) + ... + r fc (T))] (36) 



+ ^log 



and 



A fc = -| Q Pk dt+— pog(l + / fc (0))-log(l + / fc _i(0))] (37) 

where 1 + is the argument of the logarithm in (|36|) . 

Let us now compute the probability distribution function for A&. The form of (|29| ) 
implies that the ^-dependent terms in (|37|) give no contribution, since they do not contain 
the conjugate variables faj. 

We are therefore left with N — 1 Gaussian integrations over p±, . . . , p k -i, Pk+i, ■ ■ ■ , Pn 
which give the following exact result for the statistics of pk- 

Vp k exp[-—— — — / (p k (s) - X k ) 2 ds) (38) 



AD(N - 1) Jo 
where 

\ k = D(N -2k + l), k = l,...,N (39) 
The probability distribution function p(Xk', T) of the k-th Lyapunov exponent A& is then: 



1 / NT ( NT t \9 \ 

^ T > = 2 \/^(^T) eXP {-*DiN=T) ih - **> J (40) 

The Lyapunov exponents A& are statistically dependent due to the constraint (||) and their 

joint distribution function has a generalized Gaussian form []. We have thus obtained a 

complete knowledge of the statistics of the Lyapunov spectrum of the matrix P t . This has 

an essential application to the problem of the computation of the multipoint correlators of 

a passive scalar advected by a random velocity field (see the Appendix). 

3 The Gaussian distribution of the Lyapunov exponent in the N = 2 case was obtained in the context of 
the passive scalar problem in ref . B . 
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4 Conclusions 



In this work we gave a detailed exposition in the general N x N case of a functional integral 
method for the averaging of time-ordered exponentials of random matrices which has found 
several applications to the study of the statistical properties of disordered physical systems 
0, [|, [I], [| |7], H and we have shown how the statistics of the Lyapunov spectrum of a linear 
evolutionary process can be computed exactly. As a matter of fact, our method provides 
a general setting for computing the statistics of linear evolutionary systems subjected to a 
white noise force field. 

We would like to conclude with some remarks. The definition of the Lyapunov exponents 
as the logarithmic rate of growth of a k- dimensional parallelepiped (see eq. (|32J) and ref. ||) is 
the most natural from a physical point of view, e.g. in the passive scalar problem. Generally 
speaking these exponents do not coincide with the logarithms of the eigenvalues of the 
evolution matrix P t . The statistics of the eigenvalues of a similar evolution matrix was 



studied in refs. fllO|, pTT], |12i pL3l - O ur method, however, allows one to obtain a more detailed 
statistical information about the evolution of initial vectors and to compute non-trivial 
correlation functions of their components. For an application to the passive scalar problem 
see ref. ||. 

Lastly, we would like to remark that a more refined application of the functional integral 
method we described allowed to solve effectively the more difficult case of a "coloured" noise 
(see ref. @). 
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Appendix 

The method we exposed in this work has a direct application to the computation of the 
statistics of a scalar passively advected by a random velocity field. In order to illustrate this 
point we will briefly recall here the terms of the problem. For more detail see refs. || || [JJ 

mi- 

The evolution of a scalar field 0(r,t) passively advected by a velocity field v(r, t) and 
generated by a source 0(r, t) is given by 

e + v -ve = (f) (4i) 
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If we impose on (f4l|) the asymptotic condition 9(r, — oo) = we get the solution 

r+oo 

6(r,t)= <f>(R(r,t-s),t-s)ds (42) 
J o 

saying that 9(r,t) is completely determined in terms of the trajectories R(ro,i) of the fluid 
particles: 

R = v(R,t), R(r ,0) = r (43) 

Let us now take <ft and v to be random, 5-correlated in time fields . The source <fi will be 
assumed to be spatially correlated on a scale L. The velocity field v might be multi-scale, 
with smallest scale larger or of the order of L. The statistics of and v will be assumed to 
be spatially isotropic. 

Generally speaking, one is interested in computing equal-time correlators of the form 
(9(ri, O)0(r 2 , 0)) for |r 2 — r x | < L. From the isotropicity of the statistics of <j> and v it follows 
that such quantities are rotation- invariant. Moreover, (f4^ ) implies that the statistics of 9 is 
completely determined in terms of the statistics of the trajectories Q43f) . 

In order to subtract the effect or sweeping, let us chose a reference frame locally comoving 
with one of the fluid particles (see refs. [§, || 0). We can then locally linearize ([43]), obtaining 

R ~ a(t)R (44) 

where = dvj/drt, the matrix of velocity derivatives, will be taken to be a random gaussian 
process. In the general case we have a = R + S, where R is the antisymmetric part of a, 
inducing a rotation of the passive scalar blob, and S is the symmetric part, representing the 
stretching of the unit blob. We will consider the case of an incompressible fluid, so Tr S = 0. 
More specifically, let us consider the following statistics of a: 

VM [a] = Vo exp f - £ C > ( 45 ) 

/ / \ 2 

C = -(ttS 2 -^TiR 2 ] =-Tr S + iJ^R) 



2 V D r J 2 \ V D r 

Since we are interested in rotation-invariant quantities the final result shall be independent 



on D r . This arbitrariness allows one to substitute \J D s / D r — > —i, R — > iR, a — > X = S + iR 
(see ref. |I[), and thus to consider a generic traceless hermitian matrix X with averaging 
weight exp (— 2]=r / l^rX 2 ) i n the place of the generic traceless real matrix a with the 
averaging weight (551): this allows one to refer to the results of sect. 2. 

From rotational invariance it follows that the statistics of 9 is completely determined in 
terms of the statistics of the Lyapunov spectrum of the matrix P t defined by 

A = XP t , P = 1 (46) 

This reduces us to the problem studied in the preceding sections. The Gaussian statistics 
of the Lyapunov exponents agrees with an old result [17]] about the Gaussian statistics of a 
line element in a 5-correlated in time velocity field. 
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